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CALCULATION OF FORCES ON MAGNETIZED BODIES 
USING COSMIC NASTRAN 
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INTRODUCTION: 

The analogy between the equations of Heat Transfer and 

Magnetostatic# has been used by several researchers to calculate 
magnetic fields in a variety of applications and shown to give 
correct results in theoretically verifiable cases. The resultant 
fields and fluxes have been used to find forces upon moving 
charged particles in the design of particle accelerators (1) and 
upon electrical conductors in the design of electrical machinery 
(2). Magnetostatic forces, due to the action of magnetic field# 
on magnetically permeable bodies, are a less tractable problem. 
Analytical methods (3,4,5) generally applied to obtaining an 
net force on a body of iron in air have produced apparently 
conflicting expressions for force# based upon the virtual work 
method, calculation of Maxwell stresses, and use of a free 
magnetic pole model. Carpenter has reconciled the apparent 

discrepancies (5) but the expressions obtained are not entirely 
suitable for finite element analysis, as they are shown by 

Carpenter to be mathematical devices to produce an equivalent 

force on a body rather than a true representation of a force 
distribution. While in many cases an equivalent force is all that 
is required, such an approach is not appropriate to problems 
involving structural deformation. 


MAGNETOSTATIC EQUATIONS 

The equations of Magnetos t at ics are based upon forces between 
field sources. In the S.I. system of units, the Ampere is the 
current required to produce unit force between two long parallel 
conductors separated by unit distance. In the older c.g.s. e.m.u. 
system, unit force was produced between unit magnetic poles in 
vacuo at unit distance. The imaginary but convenient concept of a 
magnetic pole 1# difficult to express in the modern system of 
units . It 1# convenient to describe the properties of a magnetic 
material in terms of the applied magnetic field, H, and the 
volume magnetization, M, where the flux density, B, is given by: 

§ =yA„W-*1 Cl) 


where is the permeability of free space. This is the Kennelly 
SI system. In the Sommerfield S.I. system the magnetization is 
dimensionally the same as the field strength and is also operated 
on by the permeability constant to obtain B: 


207 



8 (2 ) 

In order to use the NASTRAN heat transfer solver the equation 
must be recast in the form : 

6 ( 5 ) 

whereyitis the permeability of the material. This formulation is 
analogous to the heat transfer equation but in no way represents 
the physical relationship between B and H. At the boundary 
between- materials of different permeability the total magnetic 
flux, v , is conserved. The normal component of magnetic flux, Bn, 
and the tangential component of the magnetic field, Hn, are also 
continuous across the interface. These conditions are essential 
to all the subsequent force derivations. 

MAGNETOMECHANICAL DEVICES 

A magnetic field has associated with it an energy given by the 
equation 

U * \ B • W 

Given that in equation (3) the value of permeability in a 
ferromagnetic material is several thousands, it is apparent that 
the energy required to support a given flux density in an air gap 
greatly exceeds that required if the same volume contains a 
ferromagnetic material. Two categories of device utilize this 
mechanically: devices employing e lec t romagnet s generate a 

mechanical force to reduce the volume of an airgap, thereby 
minimizing magnetic field energy, when a magnetizing force is 
generated by a coil, while in stored energy devices a permanent 
field is generated by a permanent magnet, producing a force that 
is eliminated by the action of an opposing el ec tr omagnet when it 
is required to remove the mechanical force. The magnetic field 
energy may be converted to kinetic energy or strain energy. 

THE VIRTUAL WORK METHOD 

The virtual work method is the most widely known means of 
force calculation. It is generally applied by calculation of 
the working gap energy, which is subsequently differentiated with 
respect to a direction of motion to obtain a force in the 
direction under consideration. In air, the usual working gap, 

8 =/u 0 H (s) 

Thus, by substitution in equation (4), an expression for air gap 
energy from which forces may be derived is obtainable in certain 
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geometries. Between two large parallel plates the force per unit 
area is: 

F = '/i/AoH 4 (fc) 


This formula neglects the energy changes within the magnetic 
material, and applies to a highly idealised orientation. Further, 
the result is a global force rather than a distributed force. 
Physically, the assumptions are that the magnetic material is of 
infinite permeability and that the source of magnetomotive force 
(MMF), is infinitely stiff. If the field energy of the magnetic 
material filling the gap is considered equation (6) becomes 

F * h/*o[\ — Jm>)U 4 CO 

ADAPTION FOR FINITE ELEMENT METHOD 

The finite element method can be applied to the virtual work 
model to obtain a higher degree of accuracy. Since the NASTRAN 
heat transfer solver calculates both B and H it is possible to 
obtain the energy associated with each element for an initial 


conf i gurat ion 
permeabi li ty 
equa t i on ( 4 ) 

and geometrical perturbations thereof 
across the element is constant, as it 
becomes simply, 

. If 
must 

the 
be , 
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If a rigid 

body is being analysed, 

perturbation in 

the 

six 


degrees of freedom, with the system energy summed in each case, 
should in theory provide complete data about the forces acting 
upon the body. The advantages of the finite element approach over 
analytical and boundary element methods are: 

(I) : Account is taken of field energy inside the magnetic materials 

(II) : The technique comprehends non-uniform fields, leakage flux 

etc. in the circuit 

Since only an equivalent force is derived the approach is 
inapplicable to problems where deformation occurs and a 
distributed force is required. It is particularly suited, 
however, to applications where the displacement or deformation is 
constrained to a single degree of freedom. Applications such as 
linear and rotary solenoidal actuators or rotating armature 
printheads may be modelled using this technique. 

THE MAXWELL STRESS METHOD 

Equation (8) may be rewritten: 
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For a spatially varying field, the equation has the more general 
form 

U s 0°) 

giving x,y and z force components of the form 

f x = ^S^ Htdxc h dft ^ 

If the permeability is constant then equation (11) has a value of 
zero. Consequently, traction forces are only exerted at the 
boundaries between regions of different permeability, or where 
the permeability varies continuously. A spatially varying field 
will not of itself produce traction forces. If a closed surface 
is constructed around a magnetized part then, from equation (11) 
a stress through the surface and a pressure tangential to it are 


obtained, 

Fs 
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The expressions 

reduce 

to zero if no discontinuities in H 

are 

enclosed. They 

may be 

resolved into components normal 

and 

tangential to 

the surface: 
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If the surface is constructed outside of a magnetic part (in free 
space) then /A .Alternatively a surface may be constructed 
within the part. In this latter case no field sources are 
enclosed, since all field sources lie at the interface between 
the part and free space, and a zero resultant force is obtained. 

ADAPTION TO THE FINITE ELEMENT METHOD 

Equations (14) and (15) may be applied by construction of a 
surface enclosing a part, located arbitrarily close to the 
boundary between the part and free space. This force acts thru a 
closed surface on enclosed field sources, and are not necessarily 
a physical representation of the true surface force distribution. 
The case will now be considered where a pill-box is constructed 
normal to the boundary, of vanishingly small thickness, and 
extending to both sides of the boundary. By this construction the 
sources of field are localized within the box and forces may be 
obtained in direction both out of and into the magnetic body: 
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where the suffixes o end i denote inside and outside values of 
the field and permeability. If yU,. is defined as the ratio of 
internal to external permeability then 

F » * i/». (i- 0?) 

Ffc » >Ao ( l >*>) H^Hto (|%) 

This appears to be a local force, and reduces to equation (7) 
for the case where the field is normal to the interface. The 
model, however, does not impart physical insight into the 
phenomena, and there are implied assumptions in the neglect ' of 
the sides of the pill box surface. 


THE FREE POLE METHOD 


The magnetization process consists of realignment of magnetic 
domains in a ferromagnetic material so as to orient then with the 
external magnetic field. Within the material, macroscopi ca 1 1 y , 
each North pole produced will be cancelled by a South pole, and 
there will be no discontinuities in H within the material. At the 
boundaries of the material, however, the alignment (polarization) 
of the domains produces free magnetic poles. The magnetization, 
M, is equal to the free pole density. The free poles are acted 
upon by the field H with a force in proportion to both H and M, 

F * H.M (W) 

The normal component of the force acting on M free poles is 
obtained by considering the force which must be applied to an 

increment dM, in order to move it to the interface from an 

arbitrary noncoincident point. Since the free poles themselves 
affect the relationship between internal and external values of 
H, an integration is required to obtain the normal force: 

F£ n - ^ a M h 

In considering the field from the poles alone, it is physically 
apparent that it must be normal to the surface on which the poles 
lie. Therefore, the tangential force is simply: 

F t = /*, H t M (21) 

These equations may be reformulated to be comparable with the 
results of the other methods discussed: 
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NASTRAN CALCULATIONS 

For the theoretically verifiable case of an air gap between two 
large parallel plates, forces on the plates have been calculated 
from NASTRAN heat transfer output by the methods described above. 

Since the methods reduce to the same theoretical base equation in 
this geometry, it is not surprising that the results are in 
agreement with theory and with each other. The problem of finding 
a non-trivial yet theoretically verifiable test problem remains. 

COMPARISON OF DERIVATIONS 

Equations for a net magnetic traction force have been obtained 
by virtual work, and Maxwell stress methods. Distributed forces 
have been obtained by the free pole method and a variation of the 
Maxwell Stress method. All methods are equivalent for cases where 
the magnetic field at an interface between the magnetic medium 
and air is normal to the surface. The equations governing the 
angle of field on either side of the medium are: 

Hi G»t<*i - HoCotMo - l*/Mo)Cd«Xc 
HC S."*;. . HoW<o 

where ^ls the incident angle of the internal field 
is the incident angle of the external field * Thus, 

=■ H o Su \ 0 / ( HoCo**b - ( ( 14 ) 

For reasonable values of M, oti* 0. Thus in all but extermely 
saturated nmagnetic materials the normal force will be orders of 
magnitude larger than the tangential forces, and the latter can 
be ignored for all practical purposes. 

Although not admired by theoreticians, the polar model of 
magnetism is mathematically valid and there is no reason to doubt 
the results obtained above. The polarization model is based on an 
analogy between elect r os t a t i c and magnetostatic phenomena which 
is not entirely valid, but which provides the same results in 
terms of force and energy as for the more correct domain model of 
ferromagnetism. The use of analogy has been compared unfavourably 
with the mathematical purity of Maxwell's equations, since it is 
little known that these were derived using a tortuous mechanical 
analog of electrostatics and magnetostatics. The Maxwell stress 
method described above is one of three stress systems presented 
by Maxwell, one of which is in total contradiction and has been 
described by Carpenter (5) as completely useless. The discrepancy 
between the Maxwell and free poles methods is likely due to the 
neglect of the side surfaces of the pill-box structure described 
above, and in any case is of negligible magnitude. 
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CONCLUSIONS 


The method* described may be used with . .. . 
confidence for calculations of magnetic traction f** h de * ree of 

• .urf.c i„ thl . eircu^tMeTSJ " 

have resulted In the theoretically correct result Tt*?!* 
abov. that the tan.snti.l force, are in J|£cti£ ? k" 
surface pole method Is preferable to the virtual Work method 

latter* case ^ U*h*"* lty f ° r "° r * th * n on * MASTR AN run in the 
latter case, and because distributed forces are obtained Th! 

derivation of local force, from th. mIJwIu iatES 

involve, an unde.ire.bl. d.*rs. of manipulation of [SI 

»thcS: oduc ** * - ,uit •» «^*«««*« », th. 
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